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Abstract 

In this paper we discuss the bihamiltonian formulation of the (rational 
XXX) Gaudin models of spin-spin interaction, generalized to the case 
of sl(r)-valued "spins". In particular, we focus on the homogeneous 
models. We find a pencil of Poisson brackets that recursively define a 
complete set of integrals of the motion, alternative to the set of integrals 
associated with the "standard" Lax representation of the Gaudin model. 
These integrals, in the case of su(2), coincide wih the Hamiltonians of 
the "bending flows" in the moduli space of polygons in Euclidean space 
introduced by Kapovich and Millson. We finally address the problem of 
separability of these flows and explicitly find separation coordinates and 
separation relations for the r = 2 case. 

AMS Classification numbers: 70H06, 37K10, 70H20. 



1 Introduction 

In [12 . M. Gaudin proved the intcgrability of iV-sitc su(2) (quantum) spin 
Hamiltonians of the form 

N 

H=Y^^^o r o u (1.1) 

where a = [a x ,a y ,a z ] are Pauli matrices. This property follows from the fact 
that one can write the Gaudin Hamiltonian 7i as 

N 



h = J2cjH 3 with ^ = E^rf ( L2 ) 
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and check that the H/s define a set of N — 1 commuting observables. Since 
fixing the values of the N Casimirs Cj = \<Ji\ 2 the system has N degrees of 
freedom, the N — 1 quantities Hj, together with, e.g., S z = ^2%=i a ti P rov ide a 
complete set of mutually commuting observables. 

A relevant member of this class of Hamiltonians is obtained when one chooses 
Cfc to be proportional to for all fc's, so that, up to a rescaling, the Hamilto- 
nian becomes 

N 

H = J2°r°i- (1-3) 

3,1=1 

This is the Hamiltonian of the XXX rational homogeneous Gaudin model. 

This system is not only integrable, but maximally superintegrable. One can 
understand this stronger property as follows (see, e.g [T5]): since the "physical" 
Hamiltonian (jl.3|) is independent of the parameters, the choice of the a^'s in the 
definition of the commuting integrals Hj is arbitrary (provided tij ^ dj, i ^ j). 
So, choosing another set of parameters bk ^ a*, and considering Hi = Tr^El 
one can define the two sets of complete commuting quantities: 

{H, Hi, . . . , -ffAr_2, S z }; {H, Hi, . . . , H N _ 2 , S x }. 

Since for generic choices of the sets a^, bk the observables 

{H, Hi, . . . , Hn-2, Hi, . . . , Hn-2-, S z , S x ] 

are algebraically independent, the model is indeed maximally superintegrable. 

Recently it was pointed out independently by various authors jU [TH] , that 
with (jl.3j) it is possible to associate a set of commuting integrals independent 
of the parameters. Such operators are of the form: 

k 

4-1 = Yl & r°i k = 2,...,N, (1.4) 

3,1=1 

and, together with, S z they form a complete set of involutive integrals for TC. 

In a completely different context the classical version of the integrals ([1.4)1 . 
namely 

^ = 2^(I>) 2 )' (1-5) 

i=i 

where A\, . . . are generic elements of su{2), were considered by Kapovich and 
Millson [T7j. These authors (see also [TT]) studied the moduli space of N + 3- 
sided polygons in M 3 , and (implicitly) showed that it coincides with a suitable 
Marsden-Weinstein quotient (with respect to the diagonal action of SU(2)) of 
the phase space of the N + 3 site sw(2)-Gaudin models. They remarked that 
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such a space possesses a natural Hamiltonian structure, and integrated, via 
action-angle variables methods, the flows associated with the integrals Ki,l = 
1, . . . , N. It is worthwhile to remind the intriguing representation of such flows: 
if one draws, from a chosen vertex, the N possible diagonals of an N + 3-sided 
polygon, the flow associated with the Hamiltonian geometrically represents 
the bending of one side of the polygon along the k-th diagonal (the other side 
being kept fixed) , whence the name of "bending flows" . 

The Gaudin system fll.lj) admits various generalizations. Gaudin himself 
pointed out that the integrals (jl.2j) can be generalized to any semisimple Lie 
algebra g. Clearly, if the rank of g is greater than 1, the number of such integrals 
is not enough to ensure complete integrability. The missing integrals have been 
shown by Jurco ^H] and Sklyanin [23] to be provided by the spectral invariants 
of a suitable Lax matrix, whose classical counterpart is 



where aj ^ dj, i ^ j and the A{ are generic elements of g. In terms of the Lax 
matrix (jl.6|) the generalization of the Hamiltonian ()1.3|) reads: 



Another straightforward generalization of this model is obtained adding a 
constant term a to the Lax matrix. In the su(2) case, this is equivalent to 
adding to the Hamiltonians (jl.3|) a term describing the interaction of the spins 
with a magnetic field with a constant direction in each site but with different 
intensity. In this case one speaks of "inhomogeneous" Gaudin magnet. The 
complete integrability and separability of these systems, (for the g = sl(n) 
case) was studied and proved in [213 EH 121] • 

The aim of this paper is to frame the analysis of the Gaudin models, as well 
as of the Hamiltonians (|1.5|) of the bending flows of Kapovich and Millson, in the 
scheme of bihamiltonian geometry as advocated by Gel'fand and Zakharevich 
[T%] , and to show how one can use this scheme to explicitly integrate the model 
for q = sl{2). We will consider only the classical models, and consider the 
complexified case (that is, we will study the Gaudin system associated with a 
complex semisimple Lie algebra g) . 

Our first task will be to briefly show how, using nowadays standard results 
of the theory of r-matrices on loop algebras (see, e.g., one can provide 

the phase space of the (inhomogeneous) Gaudin magnet with a bihamiltonian 
structure, selecting it out of a multi-parameter family of Poisson structures. 
This structure gives rise, according to the GZ scheme, to the integrals associated 
with the Lax matrix of Jurco and Sklyanin (jl.6|) . 




(1.6) 



A' 



Hg = E res U=* Tr(AL r 2 at ) = ^Tr(A^,). 



(1.7) 
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Then we will construct, in the homogeneous case, another bihamiltonian 
structure, non compatible (in a sense to be precised later) with the abovemen- 
tioned family, whose GZ analysis gives rise, in the si (2) case, to the parameter 
independent integrals (|1.5|) . Since such additional bihamiltonian structure is 
still constructed within a Lie-theoretical setting, we will be able to straightfor- 
wardly apply this scheme to g = sl(r), with arbitrary r. In this way, we will 
be able to find a sufficient number of commuting integrals to be added to the 
"generalized bending Hamiltonians" yielding a complete family of integrals 
alternative to the "standard" family obtained by Sklyanin and Jurco. 

The GZ analysis of such a model will finally lead us to introduce a kind of 
Lax matrices for such flows and to show that the Hamilton- Jacobi equations 
associated with the si (2) bending Hamiltonians are separable by computing 
explicitly the separation variables and the separation relations. 

2 GZ analysis of Gaudin models 

The Gel'fand-Zakharevich (GZ) scheme ^3] for integrating a bihamiltonian 
system can be seen as a particularly efficient scheme to implement the Lenard- 
Magri recursion for manifolds endowed with a pair of compatible Poisson brack- 
ets none of which is symplectic (i.e., non-degenerate). 

One considers a manifold M endowed with a pair of compatible Poisson 
tensors P\ — XPq, or, in other words, a pencil of Poisson brackets 

{f,g}x = {f,g} Pl - HLg} Po = (df, (Pi - \P )dg). (2.1) 

(where (■, ■) is the canonical pairing between T*M and TM), and assumes that 
the kernel of the generic element of the Poisson pencil be fc-dimensional. Let 
Ci, . . . , Ck be independent Casimir functions of Pq. The GZ method, roughly 
speaking, suggests to use these Casimirs as "starting" elements for Lenard chains 
yielding (under some technical additional conditions), via the method of biham- 
iltonian iteration, families of functions {H$}™ =1 k — 0, . . ., such that for any 
function F on M, and a — 1, ... ,k, 

{F, H$U = {F, H^} Pl , with Hi a) = C a . (2.2) 

As a consequence of the bihamiltonian iterative scheme and of the fact that 
all Lenard chains start with a Casimir function of Pq (they are "anchored" , in 
the language of [El)) all these functions are mutually in involution with re- 
spect to both Poisson brackets. Obviously, the maximal number of independent 
functions one may hope to get in this way is N max = |(dimM + k). If this 
is indeed the case, the geometric scheme herewith outlined defines families of 
completely integrable systems in the Liouville sense. Indeed, let us suppose that 
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the GZ method provides us with k families of mutually commuting independent 
functions 

{h^ }r=t::.x a , with E n - = 5 ( dimM + fc ) • 

0=1 

Let 7i be a generic element in the ring generated by such commuting functions, 
and let Xy_ = P^dTi be the corresponding Hamiltonian vector field. Let us con- 
sider a generic symplectic leaf S C M of Pq] it is a d s = dimM — k dimensional 
manifold, with the natural symplectic form induced by the Poisson structure Pq. 
X-h clearly restricts to S, and, as a consequence of the bihamiltonian iteration 
on M, comes equipped with |(dimM + k) — k — |c?5 integrals in involution, 

given by the restriction to S of the functions {Hj% Y^=\'"'k a ■ ^ s a consequence 
of the genericity assumption on the symplectic leaves, these functions will be 
independent on S as well and give the complete family of involutive integrals 
required by the Liouville theorem. 

The aim of this Section is to frame the (general, that is, inhomogeneous) 
Gaudin model within the bihamiltonian scheme, and to reinterpret its complete 
integrability within the theoretical framework of the GZ analysis briefly sketched 
above. The manifolds we will consider will be Cartesian products of a Lie algebra 
g N , and the Poisson pencils suitable linear pencils of g N . Since these results 
are essentially known in the literature, we will mostly limit ourselves to state 
results and sketch proofs, referring to [2*3*1 "2~2~| for a more general setting, and to 
[TU] for the explicit study of the 3-particle s/(2)-case. 

2.1 Notations and Conventions 

Let us briefly recall the notion of Lie-Poisson brackets associated with a Lie 
algebra and fix some notations and conventions we will use throughout the 
paper. 

If g is a Lie algebra, its dual g* comes equipped with the standard Lie-Poisson 
structure: 

{F, G}{A) = (A, [dF, dG}) = (dF, PdG), F,G G C°°(g*). (2.3) 

If g is semisimple we can identify g* with g. Indeed, we can associate a matrix 
Xa with any element A e g*. considering, e.g., the fundamental representation 
of the algebra g, and taking the trace form as a bilinear non-degenerate pairing 

(A,X B ) = Tr(X A -X B ). 

From now on we will implicitly use this identification, and write A,B... in- 
stead of Xa-iXbi ■ ■ ■ for simplicity of notation. Using the ciclicity of the trace, 
the Hamiltonian vector field associated by (J2.3j) with a smooth function F is 



5 



represented by: 



PdF = A 

where the symbol denotes the matrix satisfying, for any S in g, 



dF 
A dA 



F(A + t~) = F(A)+t.Tr(^.~) + o(t). 

If we take the direct product of N copies of g, the standard Lie-Poisson structure 
becomes: 

N r f)F fic 

{F,G}(A U ...,A N ) = J2(A, 



pAi dA-i 

and the Hamiltonian vector field associated with a function F is 

dF' 



(2.4) 



A, 



OA, 



i = l,...,N. 



We can write the above equation in the form: 



OA 
dt 



= (X F ) i = (PdF) i = ^2p ijk 



A k , 



dF 
dA~ 



with p ijk = 5 i: j5 jk . (2.5) 



We will also often write P (and other Poisson tensors) representing its action 
on the differential of a function by means of the matrix symbolic form: 



( a\ 



\A N ) 



( 



o 



[A 2 ,.\ ... 








\ 



V 



/ dF \ 

dA 1 
dF 



[A N , .} J 



dAo 



dF 



(2.6) 



V dA N ) 



For this reason, we will term the standard Lie-Poisson tensor P on g N the 
diagonal Poisson tensor. 



2.2 A Bihamiltonian Structure of the Gaudin model 

A bihamiltonian structure for rational Gaudin models can be obtained using 
the following argument. Let us consider the map {Ai} — > {Bi} that sends the 
rational Lax matrix 

v— v A, 

Wat = 0~ + 2 < 



i=l 



A - a,- 
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in the polynomial Lax matrix 



N-l 



N 



L poly = \ n a+J2 B ^ = ( I[( A - ' L 



rat 



(2.7) 



i=0 



i=i 



given explicitly by 



N 



8=1 

-l) N - l s N _i(ax,. 



SjV-i-H°i> 



, a N ) ■ a, 



,a N ) ■ Ai+ 



1 = N 



(2.8) 



where Sfc(ai, . . . , a^) denotes the fc-th elementary symmetric polynomial in the 
variables ai, . . . ,a N . 

On the space of polynomial pencils of matrices a family of mutually com- 
patible Poisson brackets are defined [231111 • They will be termed, for the sake 
of brevity, RSTS tensors. In a nutshell, this family can be described by saying 
that there is a map from degree N polynomials in the variable A to the set of 
Poisson structures on the manifold of polynomial Lax matrices of the form (|2.7J) 
which sends the monomials A , . . . , X N into N + l fundamental Poisson brackets, 
Hi, I = 0, . . . , N. In our case, the fundamental tensors II; can be represented by 
matrices having the following block-diagonal structure: 



n, 



D, 



(2.9) 



with: 



(Ci)ij = -[Bi+j-i-x, .] i,j = 1, . . . , I 

( D i)ij = [ B i+j+i-u ■} i,j = !,•••, N -I 
Bi = if % < or i > N, and B N = a. 



(2.10) 



Lemma 1 In the "coordinates" B Q , . . . , -Bjv-i, <J , the diagonal Poisson tensor 
P \2. 6)) is given by the sum 



N 



P = ^-l^-'-V-itei. • • • > o>n)Hi 



(2.11) 



1=0 



where the Sj 's are the elementary symmetric polynomials in the a, 's, that is, it 
is the tensor associated with the polynomial 



N 

Pn = JJ(A - a^. 

i=i 
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This Lemma can be proved by means of a direct computation. For the reader's 
convenience, we collect the main steps of it in Appendix A. 

Since the Poisson tensors ()2.9j) form a. N + 1-parameter family of compatible 
Poisson tensors, we can choose as a second Poisson tensor a suitable linear 
combination of them to have a bihamiltonian structure on g N . Let 

N-l 

Q = ^(-l^-'sjv-i-iCai, • • • , a N )U h (2.12) 
1=0 

be the tensor associated with the polynomial 

Pn—i = (^) + = A^ 1 - Sl A^ 2 + ■ ■ • + (-l)%v-i. 

All the integrals of motion that one can obtain from the spectral invariants of 
the Lax matrix ()1.6|) can be obtained by the GZ method applied to the pencil 
Q — \P; in fact it holds (see, also, [2H]): 

Lemma 2 All the vector fields associated to the spectral invariants of M.b}) are 
bihamiltonian with respect to the pair Q — XP. 

Proof: We find convenient to work in the variables E>i. Let us define: 

^ = Tr(Res A . (^%^)) (2.13) 
i = 1, . . . , aN a = 2, . . . rk(g). 

For any fixed a, the aN functions ()2.13|) fulfill the relations [23J: 

UidK^ = Ii i+k dK^ +k) = xU-*\ (2.14) 

From (I2~T1) it follows that X® = if i < or i > JV(a-l); in fact, if i < then 
K® = and X® = U dK® = 0, while if i > N(a - 1), then #(^+*) = const 
and X® = Il N dK( N +^ = 0. Now let us set: 



l) N - l+1 s N _ l {a 1 ,...,a N ), (2.15) 



we have: 



N N 

PdK ij) = biTlidK( j) = kX (j - l) 

1=0 1=0 
N N 

QdK® = J2 b i n i-i dKQ) = J2 blX(J ~ l ~ 1] ■ 
i=i i=i 



Then: 



PdK^-QdK^ = b Xi j l 
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If one of the aj is equal to zero, then 6 = IL=i a i = and the proof is concluded. 
Otherwise we need to find a function such that 

QdF® = boX® 

We proceed by induction. If j = 1, we have b Q X a ^ = Q^dK a , so that Fa = 
bo/biK^. Now let F® be such that: b X® = QdF$ i = 1, . . . , Then 

Q^dK? = b x^ + b -^xt 1] + ■■■ + b -^xt N+1) =► 

Oi Oi Oi 

hX® = Q ( £dK<p - b ^dFt 1] ^dFt N+1) ) ■ 

\bx b\ b\ J 

So we have: 

b 1 N ^ 

QdFW = boX®, with F® = -±K® - r J2 b i+ iFt i] (2-16) 

□ 

Some observations on the GZ sequences are in order. The starting points of the 
GZ sequences are given by the Casimirs of P. We have to distinguish two cases: 

a) If b 7^ 0, i.e. all the Oj are different from zero, then the Casimirs of P are 

given in terms of the spectral invariants (j2.13|) by the following expressions: 

aN 

C hQ = J2< Ki J ) i = l,...,N a = 2,...rk(g). (2.17) 

3=1 

For any a, starting the GZ sequences from suitable linear combinations of 
the Casimirs Cj jCe we can construct N GZ sequences of length a — 1 (i.e. 
defining a — 1 independent vector fields) each starting with a Casimir of 
P and ending with a Casimir of Q. 

b) If bo = then only one among the czj, say aN, is zero. In this case equation 

(I2.17|) defines (rk(g) — 1)(N— 1) independent Casimirs, instead of (rk(g) — 
1)N: 

aN 

= ^24 K a ] i = l,...,N-l « = 2,...rk(g). (2.18) 

3=1 

The functions ()2.18|) turns out to be simultaneous Casimirs for both P 
and Q. The remaining rk(g) — 1 Casimirs of P (the rank of P is obviously 
the same in both cases) are given by 

C Na = K^ a = 2,...,rk(g). (2.19) 

To each Casimir (J2.19j) is associated a GZ sequence of length (a — 1)N. 
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3 The homogeneous case 



The constant term a in the Lax matrix (jl.fi j) physically describes the coupling 
of the i-th spin with an "external magnetic" field (3i = at a. The matrix a in the 
definition of the rational Lax matrix (|1.6|) is somewhat a free parameter in the 
theory. Changing a amounts to "changing the direction" of this magnetic field. 
The choice usually done in the literature is the generic one (say, a is a diagonal 
matrix with different entries); this ensures the functional independence of the 
coefficients of the spectral invariants of L rat , whence the fact that they are in a 
sufficient number to yield complete integrability of the model. 

If a is not generic, but the dimension of its stabilizer g CT := {r 6 g s.t. [r, a] = 
0} is greater than the rank of g the following happens. Not all the spectral in- 
variants of the Lax matrix are functionally independent, but one can recover 
the "missing" integrals noticing that the functions: 

N 
i=l 

commute with all the spectral invariants of L rat . 

However, something more substantial occurs for a = 0, that is, in the homo- 
geneous case. As we have recalled in the Introduction, in such a case Hq (jl.7j) 
defines, for g = si (2), a superintegrable Hamiltonian system and, in particular, 
it is possible to find another complete set of commuting first integrals which are 
not explicitly dependent on the parameters a^. 

From now on we will focus on this additional family of integrals, that, in the 
classical iV-site sl{2) model can be given by the very simple formula: 

i 

h-i = Tl ( A i -A h ),l = 2,...,N. (3.1) 

j,k=i 

We will introduce a further Poisson structure R on the manifold q n . As we 
shall show it will be possible to combine it with the diagonal Poisson structure 
P of Eq. ()2.6j) to get a further Poisson pencil, not belonging to the RSTS family 
described in Section El The GZ method applied to the Poisson pencil R — \P 
will give rise to these new set of integrals. Since everything will be done in a 
Lie-algebraic setting, these results hold for a generic semisimple Lie algebra, 
and in particular, for sl(r) with arbitrary r. 
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3.1 The additional bihamiltonian pencil 



Let us consider the bivector R, defined, by means of the constructions outlined 
in Section 121 by the following matrix: 



R 



[Ai,.] 
[At,-] [A 2 -A 1 , 



[Au-\ 
[MA 



\ 



\[Ai,-] [A 2 , 



[(N-VAv-^A,-) J 



(3.2) 



Proposition 3.1 The bivector R defined by \3. is a Poisson bivector, and it 
is compatible with the diagonal Poisson tensor P. 

Proof: Linearity and antisymmetry are obvious, so we must prove only the 

Jacobi identity. Also, we can limit ourselves to prove the assertions for the 

case of linear functions on g. If F, G, H are such functions, identifying their 

dF 

differentials with the three iV-tuples of matrices {Pi}, {7i}, (e.g., 77-— = 

oAi 

«j , . . . ) , the Poisson bracket is defined by: 

{F, G} R = (dF, RdG) = r ijk Tr (a, [A k , (3,}) = £ r ijk Tr (A h [/?,-, aj) , 

i,j,k i,j,k 

fijk — (k — l)5ij5j k — 6(i_ k )5ij + 9y-i)6ik + 6(i-j)5j k , 

where 5 is the usual Kronecker symbol and 9 is the discrete Heaviside function, 
normalized as 

f 1 if % > 
%) = \0 ifi<0 

The Jacobi identity reads 

{H, {F, G} R } R + {F, {G, H} R } R + {G, {H, F} R } R = 
Y r ijk r imj (Tr (A k UPm, «/] , 7i]) + Tr (A k [k, li] , AD + Tr (A k [[7 ro , A] , Ot])) 

i ,j,k,l,m 

which, renaming the indices, becomes 

^ r ijk r lmj Ti (A k [[/3 m , a{\ ,Ji)) + r mjk r aj Tr (A k [[a h ji) , An]) + 

+ r ljk r mij Tr (A k [[y u An] , a/]) . 
A sufficient condition for the last expression to be zero is that for any k it holds: 



(3.3) 



11 



In fact, a consequence of ()3.3j) to hold is that implies that tikim = J^j r ijk r imj is 
invariant for cyclic permutations of the indices i,l,m. So, if (j3.3|) holds we can 
write: 

{H, {F, G} R } R + {F, {G, H} R } R + {G, {H, F} R } R = 

= ^2 ^™ Tr (Ah ([[^m> a l\ » Ti] + Ti] > An] + [[Tii /^m] , on))) , 

which vanishes thanks to the Jacobi identity in g. 

Let us show that (|3.3|) holds in our case. By means of algebraic manip- 
ulations, namely cycling through i,l,m and renaming the indices using the 
Kronecker's 5, we obtain: 

3 3 

= &ikblm[{l — / ^ + fiilfimk[{i ~ m )@{i-m) + / ^^(j-mj^ji-j)] + 

3 3 
+8lm[0(i-k)(d(l-i) — 0(l-k)) + 9{l-k)9(i-l)] + Sn[0(m-k)(9(i-k) — 9(i-m)) — Q{i-k)6{m-i)} + 
+Sik[0(l-i)(O(m-l) — Q(m-i)) + 0(m-i)0(i-m)] + <W|y(J-m) (%-m) _ %-J)) — %-m)0(i-i)] 

Using the identities: 

z20(j-i)0(i-j) = {l-i- l)8(i-i) 
j 

8{i-k)(8(i-i) — 6(i-k)) + 9(i-k)9(i-i) = —9(i-k)8u 
we see that every term cancels out. 

We now prove that R is compatible with the diagonal tensor P. The character- 
istic condition for the compatibility of two Poisson tensors is 

{H, {F, G} P } R + {H, {F, G} R } P + cyclic permutations ofF, G, H = 0. 



dG dF 
dA, ' 8Ai 



one shows that a suf- 



Recalling that {F, G} P = J2i,j,k 4/<^fc Tr M.* 

ficient condition for the compatibility of R and P is that the quantity 

Siklm ^ (j ijk&lm&mj ^ij^jk^linj) 
j 

be invariant under cyclic permutations of the indices i, I, m for all k's. Actually, 

Siklm — (k + i)SikSi m 5u — 0(i_ k )5i m 5u + 9^ m )Si m Skl + 0(m-l)^ik^il + 9(l- m )^ik^im 

that manifestly satisfies this property. 

□ 
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Remark. By the previous proposition, we can endow, for every N, the space 
(q*) n with a bihamiltonian structure P\ = R — \P. A natural question arises, 
that is what is the connection with the RSTS family of Poisson structures 
discussed in Section El We do not have yet a complete answer to this point; 
however, as we will show in Appendix B, the new tensor R is not compatible 
with the generic element of the RSTS family ()2.9|) . 



3.2 The Lenard Chains 

We construct the Lenard chains for the Poisson pencil R — XP, using the GZ 
recipe discussed in Section El To shorten notations we define: 

i-i 1 
Bi = Y, **> F i? = reSA=0 ^+T ( AA < + B if ( 3 - 4 ) 

i=l 

K) = Tr (^?)- ( 3 - 5 ) 

We will first find a kind of "modified" recursion relation. 
Lemma 3 It holds: 

PdHf~ l) = (R-(l- 2)P)dH%> (3.6) 

fora=l,...,P 

Proof: We proceed by induction on /3. If /3 — 2 we have: 

H$> = 2tr(AtBi) =► ^ = 2fy_i)4 + Bfa 

H$ = tr(Af) =► ^ = 2^1. 

By direct computation we obtain: 

{{R -(I- 2)P)dH®)i = 2(d {l _ t) [A t , A t ] + duiA, B,}) = (PdH^)i 
((R -(I- 2)P)dH^% = 20 {M [A, Bi] = {PdH^% 
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We use the inductive hypothesis in the case a < [3 — 1. Plugging in the following 
identities: 

[B l ,F^X l ] + [A l ,F^ l > ] = (3.7) 
' Vi)^ + ^ ( 3 - 8 ) 



dAj K J ' dBi 1 dA 



if = "^T (3 ' 10) 

and using the inductive hypothesis, one obtains by straightforward computation 

~(R-(l- *)P)dH$ = 0. 

The case a = (3 can be easily verified by direct computation. 

Notice that identities ()3.7|) and ()3.8|) follow from (|3.4|) and (|3.5|) . The identity 
(|3.9|) follows from the recursive formula for Fpf: 

*PJL ~ °l b (3-1,1 + A l*i3-l,l 

= if a > /3 or a< 0, 
while (|3.10|) can be proven again by induction. For a = 2 it holds: 

<9# (2) «9# (1) dH (l) dH {0) 

— 2,; u 2,1 _ 2B uri 2,i 



dA t dB t dAi dBi 

Then, if a < (3 — 1 we can use the inductive hypothesis and get: 

dBi " l [ dAi dBi ) + l \ dAi dBi 1 " 

The case a = (3 is again a matter of simple computation. 



□ 



Proposition 3.2 The Hamiltonians 



k - 1 'k - 1 



3=0 



J 



Kf} = H)S, 1 = 2,. ..,N 
satisfy the standard Lenard-Magri relations PdK^ = RdK^l, a = 1, ... (3, V/3. 
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Proof: Using Lemma we have: 



RdK 



09-*) 

13,1 
k-l 



k-l 



((I - 2) k -^ 1 Hf- j ~ 2) + (/ - 2) k -^Hf- J - l) ) 



Pd K 



((3-k-l) 
(3,1 



k-l r 

E 



k-l 



fc- 1 

J 



(I - 2) k - ] H 



k-jfjiP-j-l) 

/3,l 



09-fc-l) 



□ 



3.3 Complete integrability for = sl(r) 

We now prove that the Hamiltonians (|3.5|) together with additional integrals 
one can recover from the global SL{r) invariance of the model provide complete 
integrability in the case g = sl(r). We start by remarking that the content of 
Lemma El and Proposition 13.21 can be rephrased as follow: If we introduce the 
N matrices: 



Li = A u 



A A* + B a , 



N 



(3.11) 



then they evolve isospectrally along any of the vector field of the hierarchy, that 
is, (since the matrices are generically simple) along Lax type equations. 

The dimension of the manifold M = sl(r) N is du = (^ 2 — 1)^, an d the 
dimension of the generic symplectic leaf of P is d$ = dM — N(r — l) = r(r — l)N. 
We notice that the we recover (as expected) all the Casimirs of P considering: 

a) the spectral invariants of L\ = Ai (this gives N — 1 common Casimirs), and 

b) the higher order terms in expansions of L\. H^, a = 2, . . . , r, I — 2, . . . , N. 



bmce it holds H$J 



E 



i,^Hpl_ l , we consider the set 



r a 



This provides us with a distinguished sequence of 



r(r— 1) 



mutually commuting 



Hamiltonians. Clearly, if / ^ the sets {H^} and {H^,} are functionally 
independent, since they depend on a different set of variables. So, the counting 
of the number of independent Hamiltonians boils down to compute the counting 
of independent coefficients in the determinant 



det(/i - \A + B), A, B e sl{r). 
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This problem was solved in j^j, (see, also, [T]) and, actually, the number is 
exactly r ( r ~ 1 > . Hence, the Lenard sequences associated with R — \P provide us 
with a total number of (iV — 1) r ( r ~ 1 ' commuting Hamiltonian, plus the iV(r — 1) 
Casimirs. For complete integrability we are missing r(r — 1)/2 more commuting 
integrals. 

They are associated with the global SL(r) invariance of the problem, and, 
in the bihamiltonian picture, can be described as follows. For every r e sl(r) 
we can consider the linear function 

N 

H T = tT(J2^ T )- 

1=1 

The Lenard "sequence" associated with such functions is somewhat peculiar; 
indeed, since RdH T = (N — l)PdH T , we can associate to each H T a Lenard 
diagram which is (up to a constant) a closed loop, to be compared with the 
usual ladder typical of iterable Hamiltonians. However, the usual argument of 
bihamiltonian recurrence, shows that, for any r, 

{K$,H T } P = {Kfl H T } R = 0, V a, (3, 1. 

Indeed, one has, e.g., the equality: 

{K$, H T } P = H T } R = (N - 1) ■ {K^\ H T } P . 

This argument shows how to recover, in the bihamiltonian formalism, the in- 
tegrals associated with the global SL(r) invariance of the model. Clearly, this 
family of r 2 — 1 integrals is not a commutative one. 

To recover the maximal Abelian subalgebra inside the Poisson algebra gen- 
erated by the functions H T , one can consider (see, e.g., [5]): 

a) The r — 1 independent elements H^, . . . , i?h r _i associated with, say, the 

standard Cartan subalgebra of sl(r); 

b) The Gel'fand-Cetlin invariants, that is, the Casimirs of the nested subalge- 

bras 

sl(2) C sZ(3) C ••• C sl(r), (3.12) 

under the map sl(r) N — * sl(r) sending the iV-tuple {A\, . . . , A^} into the 
total sum, A tot = Y^iLi 

Noticing that the Gel'fand-Cetlin functions corresponding to the last element of 
the chain f)3.12|) are given by Ylk=o HpNi we obtain r — 1 + YJiZo, — -0 = 
additional commuting integrals, which is exactly the number of commuting 
integrals we were looking for to insure complete integrability of the model. 

We end this Section with a comment concerning super-integrability of the 
model. To this end we remark that we have at our disposal two Poisson pencils 
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to construct families of commuting integrals for the Gaudin (homogeneous) 
Hamiltonian Hq\ the pencil R — XP and the pencil Q — XP, described in 
Section I2] 1 . On the d>N,r = N(r(r — l))-dimensional generic symplectic leaves of 
P they give rise to two distinct djv, r /2) families of integrals of the motion K 1 ^ 
and K%. Direct computations (which we performed for r = 3,4 and N < 6) 
suggest that the number of functionally independent elements in the union of 
the two families be d,N, r — [r — 1)- In other words, also taking into account 
the integrals coming from the global SL(r) invariance of the model, we have 
super-integrability for the sl(r) Gaudin model, that, however, is maximal only 
for the sl{2) case. 



4 Separation of Variables for the sl(2) case 



We consider now the iV-particle si (2) case. The aim is to show that the 
Hamilton-Jacobi equations associated with the Hamiltonians 



a-l 



H a = Ti(A t -J2^)- 



(4.1) 



and, in particular, the H-J equations associated with the physical Hamiltonian 
Hq = | X^j=i Hi are separable in a very "simple" set of coordinates. Our analy- 
sis is based on the so-called bihamiltonian scheme for SoV, recently introduced 
in the literature (see, e.g., [2HI2IZ1)- in particular, we will use the results for 
systems with an arbitrary number of (anchored) Lenard chain exposed in 0. 

We consider the manifold M = sl(2) N , endowed with the Poisson pencil 
R — XP, explicitly parametrized with the iV matrices 



A, 



hi fi 
Ci hi 



(4.2) 



The generic symplectic leaf S of P is a 2N dimensional symplectic manifold, 
defined by the equations 



-TrA? = h 2 i+ ei/, 



1 N, 



and can be (generically) endowed with the 2N coordinates {hi, /i}i=i,...,jv. 

According to the bihamiltonian scheme jHj, one modifies the tensor R in 
order to obtain a second compatible structure on S. Let us define the N vector 
fields 

z - ^ j,h (43) 

1 With the proviso in mind that one has to set a = in those formulas 
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and, for any pair of functions F, G on M, the following brackets: 



N 



{F,G} R , = {F,G} R -Y1 {{F,H a } P Z a {G) - {G,H a } P Z a (F)), (4.4) 



a=2 



where X(F) denotes the Lie derivative of F with respect to the vector field X. 

Thanks to the specific form of the vector fields Z a , these new brackets restrict 
to S and give rise to Poisson brackets on S, which are compatible with the 
ones naturally induced by P. So, S is a symplectic manifold with respect to 
the restriction of the brackets associated with P, and is endowed with a (1, 1) 
tensor TV, with vanishing Nijenhuis torsion defined by 

M = R'oP- 1 . 

In such a geometrical setting, the bihamiltonian scheme for SoV considers sets 
of coordinates {ui,Vi} (called Nijenhuis coordinates) associated with the eigen- 
values Aj of Af, characterized by the equations: 

M*dui = Xidui, Af*dvi = Xidvi, (4.5) 

whose Poisson brackets attain the remarkable form |2T)| ITlj : 



{Ui,Uj} P = {Ui,Uj} R i = {Vi,Vj}p = {Vi,Vj} R/ = 

{ui, Vj} P = Sijtfiiui, Vi), {ui, Vj} R > = Xi{ui, Vj} P , 



(4.6) 



for some functions Vi). 

We shall prove the above statements directly displaying a set of Nijenhuis 
coordinates on the symplectic leaves of P. 

Proposition 4.1 The 2N functions 

Ai = ^ fi, \ a = - ^ k f Jk +(a-2), a = 2,...,N 
i—i J a 

(4.7) 

Hx = ^^hi, fx a = (A a - (a - 2)) h a + ^h k , a = 2,...,N 

i=l k=l 

provide a set of Nijenhuis coordinates on S . In particular the A a , a = 2, . . . , N 
are the non-vanishing eigenvalues ofAf*, while Ai and \i\ span its (2- dimensional) 
kernel. 

Proof. According to Equation (|4.4j) and the definition of the Nijenhuis tensor 
A/", noticing that both P and R' restrict to S and that Z a (Xi) = Z a (fii) = for 
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a = 2, . . . , N, % = 1, . . . , N, we have to show that, for any coordinate Xi, (that 
is, Xi = e h h h f h i = 1, . . . , N) it holds 



N 



{xi, Xi}r - ^{H a , \i}pZ a (x,j) = 



a=2 
N 



(4.t 



{^i, A fc }i? - y^Hg, \ b } P Z a (xi) - \ b {xi, \ b } P , b = 2,...,N 



as well as 



a=2 



N 



{x u Hi} R - ^{H a , fi 1 } P Z a (x i ) = 



a=2 

N 



(4.9) 



{x^ ii h } R - ^{H a , fi b } P Z a (xi) = X b {xi, fx b } P , b = 2, . . . , N. 



a=2 



The proof that these equations holds true is a matter of direct computations. 
One simply has to plug the explicit expressions of the Poisson brackets 



{hi,ej} R = 5ij 
{hi, e.j}n = 
{ei, fj}R = 2< 8ij 



{hi, Cj}p — Sijtj, {hi, fj}p — —Syfj, {ei, fjjp — 2Sijhj 

i-l 

l)e* - y^e fc 

k=l 

i-l 

k=l 

i-l 

i - l)hi - y~]hk 
fc=i 



+ 0(i-j)ej + 0{j-{)ei 

°(i-j)fj ~ 6 {j-i)fi 

+ 0(i-j)hj + 9(j~i)hi 



into equations ()4.8|) and f)4.9j) . and use the identities 



n-l 



n-l 



n-l 



n-l 



n-l 



i-l 



j=i j=i j=i 



(4.10) 
(4.11) 



For example, let us consider Xi = hi. Since Z a {hi) = 0, we have (for n > 2): 

fi 

{hi, X n } P = (A n — n + 2)5 in + 9(n-i)-r, 

In 
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and 



n-1 



{hi, \ n }n = -J- ' S ^ J \ & 
U .1=1 



i-1 



k=l 



+ 0(i-j)fj + > + 



A n - n + 2 

fn 



i-1 



fc=l 



+ @{i-n)fn + 9(n-i)fi 



A„ Ui„(A„ - n + 2) + 0(„_i) 



thanks to identities (jOHf . (|Q!T|l . 

The other cases of equations (|4.8|4.fjjl are proved with similar computations. 

□ 

To construct a set of canonical Nijenhuis coordinates (usually considered in the 
bihamiltonian scheme for SoV and quite naturally termed Darboux-Nijenhuis 
coordinates) 4>i}i=i,...,N from the Nijenhuis coordinates {Aj, /ii}i=i,...,Ar one 
notices that an explicit computation gives: 

{A 1 ,/i 1 } = -A 1 , {A a ,/i a } = (A a -(a-2))(A a -(a-l)),a = 2,...,iV. (4.12) 

Indeed, the first equation is trivially verified; for the remaining set of iV — 1 
relations one has 



72—1 n— 1 



{fin, A„} = — T-{ K ^2 fi} ~ ( A n 
Jn i=l i=l 
Y^n— 1 n \r~\n— 1 /■ 

_ 2^=i /< _ (An _( n _ 2) )Ai=i /< . 



n-1 



in 



2)){h n ,l/f n }J2fi 



i=l 



fn 

Hence, one can choose 



/n 



V x ' 



(A a -(o-2))(A a -(a-l)) 



(A n -(n-2))(A n -(n-l)). 



a = 2,...,JV (4.13) 



to have, together with the A, a set of Darboux-Nijenhuis coordinates. 

To find the separation relations, we make contact with the so-called Sklyanin 
magic recipe [27]. To this end, we modify the Lax matrices (|3.11|) . by a suitable 
shift in the spectral parameter A. Namely we define 



a-1 



L X = L X - L a =(\-(a-2))A a + J2 A b, a=2,...,N. 



(4.14) 



6=1 



We notice that the spectral invariants of L a are combinations of the Hamil- 
tonians H^f and of the common Casimirs we considered in Section 13.31 and 
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provide an equivalent set of involutive constants of the motion (together with 

As one can easily notice, the Nijenhuis coordinates A a , a = 2, . . . , N are 
nothing but the zeroes of the (1,2) entry of the Lax matrix L a of eq. (|4.14jl . 
while the Nijenhuis coordinates fM a ,a = 2, . . . , N are the values for A = A a of 
the (1,1) entry of the same matrix. Taking into account that \i\ is the first 
Hamiltonian Hi, we see that the Nijenhuis coordinates, the Hamiltonians and 
the Casimirs C, = Tr(Aj) 2 are related by the separated equations: 

o-l 

tMi-Hi = 0; Det(/i a - ((A - (a - 2)L4 a + ^ A)) = 0, (4.15) 

i=i 

whence one can directly find, using the relations (|4.13jl . the "canonical" sep- 
aration relations for the Hamilton- Jacobi equation associated with the si (2) 
Gaudin Hamiltonians. We notice that these are quadratic equations in the sep- 
arated coordinates, and hence explicitly solvable by elementary functions, for 
every number of sites N. 

5 Summary of the Results 

In this paper we have used tools from bihamiltonian geometry to study the 
integrability of (rational XXX) Gaudin models, associated with the Lie algebra 
sl(r). We first framed the general (that is, inhomogeneous) model within the 
Gel'fand-Zakharevich scheme by selecting a suitable pencil of Poisson brackets 
induced by a natural family of Poisson brackets on the space of matrix polyno- 
mials. 

Then we extensively studied the homogeneous case. We considered an alter- 
native complete set of mutually commuting constants of the motion Ik (inde- 
pendent of the parameters usually entering the formulation of Gaudin models). 
These integrals (actually, a subfamily thereof), in the su(2) case, coincide with 
the Hamiltonians of the bending flows of Kapovich and Millson on the moduli 
space of polygons in the Euclidean space. 

We introduced an "additional" Poisson tensor which forms, together with 
the standard Lie-Poisson tensor, a bihamiltonian pencil. The GZ analysis of 
such a bihamiltonian structure provides exactly the alternative set of constants 
of the motion Ik for the sl(2) case. By using such a bihamiltonian scheme, 
we extended this analysis to the sl(r) case; in particular, we show that the 
higher rank counterparts of the additional set of integrals guarantee complete 
integrability of the s/(r)-Gaudin magnet. Since, we still have at our disposal 
the Jurco-Sklyanin integrals, we conclude that the Gaudin magnet is super- 
integrable (although we could not establish maximal super-integrability) also in 
the sl(r) case. 
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We furthermore have explicitly shown in the si (2) case that the Hamilton- 
Jacobi equations associated with the set of additional integrals can be solved 
by separation of variables, using the bihamiltonian scheme for SoV which has 
recently been considered in the literature. Actually, what we found is a set 
of separation coordinates alternative to the "standard" one found by Sklyanin 
and the "Montreal group", based on the standard Lax representation for the 
(homogeneous) Gaudin model. This should not be regarded as a surprise, since 
it is well known that superintegrability is related with the existence of different 
sets of separation coordinates. In this set of coordinates, the H-J equations can 
be explicitly solved by elementary functions, or, otherwise stated, the separation 
coordinates (for the s/(2)-magnet) live on genus spectral curves for any number 
of particles, while, in the "standard picture", the genus of the spectral curve 
grows linearly with N. 
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Appendix A 

In this appendix we sketch the proof of Lemma 1. Actually we will prove the 
converse statement, i.e., that under the map ()2.8|) the Poisson tensor ()2.11|) is 
sent exactly in the diagonal Poisson tensor (|2.6|) . 

We denote with J the Jacobian of the transformation: 

Jij = ^q^ 1 = (-l) N ~ ls N-i(ai, a N ). (A.l) 

Using the identity: 

N 

y^V"^-!)^ s N -j(ai, ...,d k ,...,a N ) = JJ(x - a k ) (A. 2) 

the inverse matrix of IA. II is easily obtained: 

(J_1) - = TT la a V (A ' 3) 
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We have: 

(J^PJ- 1 )*), 

with: 



:-i) 



N 



r in in / 



(A.4) 



N-l N 



p in =EE (-l)W-i(ai,...,a JV ) < +l - k <~ l [B r ,.\ (A.5) 



r=0 Z=r+1 
JV r-1 



fc=r+l 
r 



^ = EE(- 1 M^-' fl *) E a^-*^- 1 ^,-]- (A.6) 

r=l i=0 fc=i+l 

Subtracting (|A.5|) and (|A.6J) . by using induction and the identities 

Sj(ai, . . . , ajv+i) = Sj(ai, . . . , ajv) + a/v+iSi-i(ai, . . . , a N ) (A. 7) 
Si(ai, . . . , a N ) = 0ifi<0ori>A 



JV 



JV 



5^(-l) , s iV _i(a 1 , . . • , a N )x l = (-1) N Y[(x - Oi) 



(A.i 



1=0 



i=l 



one proves that the coefficient of B r in formula (|A.4|I vanishes if % ^ n. 
Now let us consider the diagonal terms. 



n / N 



Pi? ~ 4 2) = E ( < £(-l)W«i, • • • , - r)^- 1 ^ 



r=0 \ 1=0 



Since 

JV 



d 



£(-1) W*K - - - , a*) J a^ 1 = (-1)"^ ( II 



i=o 



we get (using (I2~%1) and (lA~2lO : 



JV 



p£ ] - p^ = (-v n - °i) e ^[Br, ■} = (-if ( - °i) ) ^> •] 



r=0 



whence the assertion. 



□ 
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Appendix B 



Lemma 4 The Poisson tensors Q ( PQ|| ) and R hS. gj) are not compatible for 
N > 3 for any choice of the parameters a\, . . . , ajv- 



Proof: An explicit computation shows that the Poisson tensor Q can be written 
in the "coordinates" {A\, . . . , An} as: 

<9G <9F 



{F,Gf} Q = ^2 qijk Tr[A k 



with: 



(1 - 5j fe ) \ (A^-jfe - ^-^fe) 



Vjk 




+ 



k=ii 



Vik 



ft 

Vji 



Using this expression is easy to evaluate the Schouten bracket of Q, R on the 
differentials of the functions 

F = Tr(Ai/i) G = Tr(A 2 x) H = Tr(A 2 h) 

where with x and h we denoted two constant matrices satisfying [h,x] = x. We 
have: 

[Q, R] s (dF, dG, dH) = 



-1) 



N 



(a _ a \ N 1 N 1 

6 - 6 + lPl P2J * r (^i*) — *r(^)(B.l 



A necessary condition for ()B.1|) to vanish is that (3\ = 0, i.e., one of the constants 
a 2 , . . . , O/v must be equal to zero. Let us suppose a fc = 0, k > 2; then 

[Q,R]s(dF,dG,dH) = {-1) N (— - — ) Tr^x) ^ 

\a 2 ai/ 

since ai 7^ a 2 . 

In the case a 2 = 0, we have instead: 

N 1 

[Q, R]s{dF, dG, dH) = (-l) N Tr(Axx) ^ — • (B.2) 



a,- 
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If iV = 3 then (|B.2|) is different from zero, so we can assume N > 3. But if 
a 2 = and iV > 3 we can consider 

F' = Tr(Ax/i) G' = Tr(A 2 x) H' = Tr(A 3 h). 

Since 

[Q,R]s(dF', dG',dH') = (-l)^ 1 — tr{A lX ) ^ 

^3 

the proof is concluded. 

References 

[1] M. Adams, J, Harnad, J. Hurtubise, Darboux Coordinates and Liouville- 
Arnold Integration in Loop Algebras, Comm. Math. Phys. 155 (385-413), 
(1993) 

M. Adams, J, Harnad, J. Hurtubise, Darboux Coordinates on Coadjoint 
orbits of Lie Algebras, Lett. Math. Phys. 40 (41-57), (1996). 

[2] M. Blaszak, On separability of bi-Hamiltonian chain with degenerated 
Poisson structures. J. Math. Phys. 39, 3213-3235, (1998) 

[3] A. Ballesteros, O. Ragnisco, Classical Hamiltonian systems with sl(2) coal- 
gebra symmetry and their integrable deformations, J. Math. Phys. 43, 
954-969, (2002) 

[4] A. Ballesteros, M. Corsetti, O. Ragnisco, N-dimensional classical integrable 
systems from Hop f Algebras, Czech. J. Phys. 46, 1153-1165 (1996). 

[5] Yu.A. Drozd, S.A. Ovsienko, V.M. Futorny, On Gel' fand-Zetlin modules, 
Proceedings of the Winter School on Geometry and Physics (Srni, 1990), 
Rend. Circ. Mat. Palermo (2) Suppl. No. 26, 143-147, (1991). 

[6] Dubrovin B.A. Matrix finite-zone operators, J. Sov. Math. 28, 20-50 
(1985). 

[7] G. Falqui, F. Magri, M. Pedroni, Bihamiltonian geometry and separation 
of variables for Toda lattices, J Nonlinear Math. Phys. 8 Suppl., 118-127 
(2001). 

[8] G. Falqui, M. Pedroni, On a Poisson Reduction for Gel'fand-Zakharevich 
Manifolds, ( |nlin.SI/0204050l ), Rep. of Math. Phys., 50, 395-407 (2002). 

[9] G. Falqui, M. Pedroni, Separation of variables for bi-Hamiltonian systems, 
|nlin-si/0204050| to appear in Math. Phys. Anal. Geom (2003). 



25 



G. Falqui, F. Musso, Bi-Hamiltonian Geometries of Gaudin Models: a Case 
Study, Proceedings of the International Conference SPT 2002 Symmetry 
and Perturbation Theory, World Scientific, Singapore, 42-50 (2002). 

H. Flaschka, J. Millson, The moduli space of weighted configurations on 
projective space, math.SG/0108191 



M. Gaudin, Diagonalisation d' une classe d' hamiltoniens de spin, J. de 
Physique 37, 1087-1098 (1976). 

See, also, M. Gaudin, La Fonction d' Onde de Bethe, Masson, Paris (1983). 

M.I. Gekhtman, Separation of variables in the classical SL(iV) magnetic 
chain, Comm. Math. Phys. 167 (3), 593-605, (1995). 

I. M. Gel'fand, I. Zakharevich, On the local geometry of a bi-Hamiltonian 
structure. In: The Gel'fand Mathematical Seminars 1990-1992 (L. Corwin 
et al., eds.), Birkhauser, Boston, 51-112 (1993). 

I. M. Gel'fand, I. Zakharevich, Webs, Lenard schemes, and the local geom- 
etry of bihamiltonian Toda and Lax structures Selecta Mathematica, new 
series, 6 131-183 (2000). 

J. Harnad, O. Yermolaeva, Superintegrability, Lax matrices and separation 
of variables, |nlin.SI /0303009| 

B. Jurco, Classical Yang-Baxter equations and quantum integrable systems, 
J. Math. Phys. 30, 1289-1293, (1989) 

M. Kapovich, J. Millson, The symplectic geometry of polygons in Euclidean 
space,3. Differ. Geom. 44, 479-513 (1996) 

V. Karimipour, Integrable Structure of the New Calogero Models, J. Math. 
Phys. 39, 913-920, (1998). 

F. Magri, C. Morosi, Quaderno S 19/1984 del Dip. di Matematica 
deWUniversita di Milano, (1984). 

F. Magri, Geometry and Soliton Equations, Atti. Acc. Sci. Torino, 124 
(Suppl.), 181-209 (1990). 

C. Morosi, G Tondo, Quasi-Bi-Hamiltonian systems and separability. J. 
Phys. A: Math. Gen. 30 (1997), 2799-2806. 

M. Pedroni, P. Vanhaecke, A Lie algebraic generalization of the Mumford 
system, its symmetries and its multi-Hamiltonian structure, Regul. Chaotic 
Dyn. 3, 132-160, (1998). 



26 



[23] A.G. Reyman, M.A. Semenov-Tian-Shansky, Group-Theoretical Methods 
in the Theory of Finite-Dimensional Integrable Systems in Dynamical Sys- 
tems VII, Springer (1994). 



[24] D.R.D. Scott, Classical functional Bethe ansatz for SL(iV): separation of 
variables for the magnetic chain, J. Math. Phys. 35, 5831-5843, (1994) 

[25] E.K. Sklyanin Separation of Variables in the Classical Integrable SL(3) 
Magnetic Chain, Comm. Math. Phys. 150, 181-192 (1992) 

[26] E. K. Sklyanin, Separation of Variables in the Gaudin Model, J. Sov. Math. 
47:2, 2473-2488 (1989). 

[27] E. K. Sklyanin, Separation of Variables: new trends Progr. Theor. Phys. 
Suppl. 118 35-60 (1995). 



27 



